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for  the  first  time.  Measurements  indicate  the  trapped-electron  belt  to  be  localized 
in  radius  and  have  a  non-Maxwell ian  energy  distribution  ranging  from  10  to  40  keV. 

Using  multiple  probes,  we  have  clearly  identified  drift-resonant  instabilities 
leading  to  rapid  radial  transport.  Transport  in  a  dipole  appears  to  require  multiple 
modes,  and  its  "bucs^y"  nature  suggests  a  profile  relaxation  of  the  energetic  electrons 
which  self-stabilLt^:^S^  the  drift-resonant  instabilities.  Theory:  Substorms  in  the 
magnetosphere  cause  the  generation  of  major  electromagnetic  disturbances  and 
energetic  particles.  We  examine  the  role  of  the  collisionless  tearing  instability  as 
a  possible  mechanism  for  substorms.  Global  asymptotic  magnetotail  equilbria  which 
are  slowly  varying  in  the  Earth-Sun  direction  are  constructed,  including  all  three 
components  of  the  magnetic  field.  Some  of  these  equilibria  are  analyzed  for  stability 
with  respect  to  collisionless  electron  tearing  modes.  It  is  found  that  the  ion  tearing 
instability,  which  has  been  widely  invoked  as  a  possible  trigger  for  subtorms,  does 
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SUMMARY 

Experiment  :  An  energetic  electron  belt  has  been  created  in  a  laboratory  terrella  for  the 
first  time.  Measurements  indicate  the  trapped-electron  belt  to  be  localized  in  radius  and 
have  a  non-Maxwellian  energy  distribution  ranging  from  10  to  40  keV.  Using  multiple 
probes,  we  have  clearly  identified  drift-resonant  instabilities  leading  to  rapid  radial 
transport.  Transport  in  a  dipole  appears  to  require  multiple  modes,  and  its  “bursty” 
nature  suggests  a  profile  relaxation  of  the  energetic  electrons  which  self-stabilizes  the 
drift-resonant  instabilities. 

Theory  :  Substorms  in  the  magnetosphere  cause  the  generation  of  major  electromagnetic 
disturbances  and  energetic  particles.  We  examine  the  role  of  the  collisionless  tearing 
instability  as  a  possible  mechanism  for  substorms.  Global  asymptotic  magnetotail 
equilbria  which  are  slow'ly  varying  in  the  Eanh-Sun  direction  are  constructed,  including 
all  three  components  of  the  magnetic  field.  Some  of  these  equilibria  are  analyzed  for 
stability  with  respect  to  collisionless  electron  tearing  modes.  It  is  found  that  the  ion 
tearing  instability,  which  has  been  widely  invoked  as  a  possible  trigger  for  subtorms,  does 
not  exist.  The  By  field  is  demonstrated  to  have  a  destabilizing  effect  on  electron  tearing 
modes.  Regimes  in  which  collisionless  tearing  modes  can  grow  are  delineated. 


THE  COLLISIONLESS  TERRELLA  EXPERIMENT 


The  “Collisionless  Temella  Experiment”,  or  CTX,  is  a  relatively  new  laboratory 
experiment  built  at  Columbia  University  in  order  to  directly  observe  the  unique  and 
fundamental  properties  of  collisionless  radial  transpon  of  plasma  trapped  within  planetaiy 
magnetospheres.  Collisionless  radial  transport  in  a  dipole-confined  plasma  occurs  only 
when  wave-particle  interactions  are  sufficiently  intense  and  broad-banded  and  when  the 
particle  motion  satisfies  well-defined  conditions  for  chaos.  The  primary  goals  of  our 
experiments  are  (1)  to  directly  observe  the  conditions  required  for  the  onset  of  radial 
transport,  (2)  to  study  the  evolution  of  the  plasma  profiles  undergoing  radial  transport, 
and  (3)  to  develop  and  test  a  generalized  model  for  radial  transport  applicable  to  trapped 
plasma  in  the  earth’s  magnetosphere.  This  laboratory  program  is  particularly  exciting 
since  (1)  CTX  is  the  first  collisionless  dipole  laboratory  experiment  capable  of  observing 
collisionless  radial  transport,  (2)  plasma  waves  can  be  produced  in  the  laboratory  which 
break  only  a  particle’s  third  adiabatic  invariant  making  radial  transport  relatively  simple 
to  characterize,  and  (3)  the  impact  of  collisionless  transport  on  a  plasma’s  global  profile 
is  predicted  to  have  an  unique  .signature  that  should  be  readily  identified. 

Our  first  experiments  with  CTX  have  focused  on  the  dynamics  of  an  “artificial 
radiation  belt”  consisting  of  trapped  electrons  with  energies  between  10  and  30  keV.  The 
electrons  are  created  by  adjusting  our  microwave  plasma  source  for  direct  cyclotron 
heating  of  magnetically-trapped  electrons.  The  electrons  generated  have  a  collisional 
mean-free-path  longer  than  1,000  drift  orbits  about  the  equator,  and  they  provide  the  ideal 
laboratory  medium  with  which  to  study  collli.sionless  radial  transport  in  dipole  geometry. 

Figure  1  shows  the  time  history  of  a  typical  CTX  discharge  and  the  production  of  an 
energetic  electron  belt.  The  discharge  duration  was  arbitrarily  programmed  to  last 
approximately  0.5  sec,  and  the  several  signals  represent  the  time  history  of  the  ion- 
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saturation  current  to  a  Langmuir  probe,  the  hydrogen  gas  pressure,  the  forward  and 
reflected  microwave  power,  and  the  x-ray  spectra  as  recorded  by  a  krypton  proportional 
counter.  The  experiment  is  fully  computer  controlled,  and  the  gas  pressure,  pulse  length, 
and  heating  power  can  be  programmed  independently  for  each  discharge.  For  discharges 
similar  to  that  shown  in  Figure  1,  intense  fluctuations  are  observed  both  during  the 
microwave  heating  and  during  the  “afterglow”  when  the  heating  power  has  been 
switched-off.  These  fluctuations  only  occur  in  the  presence  of  the  energetic  electron  belt, 
and  we  have  characterized  them  as  a  complex  and  nonlinear  development  of  drift- 
resonant  hot  electron  interchange  instabilities  (HEI).  In  some  respects,  the  HEl  is  the 
electrostatic  “analog”  of  symmetric  electromagnetic  kinetic  Alfven  instabilities  occurring 
in  planetary  magnetospheres. 

Figure  2  shows  the  probe  fluctuations  on  a  faster,  10  msec,  time-scale.  During  the 
microwave  heating,  quasi-periodic  “bursts”  are  observed:  whereas,  during  the  afterglow, 
the  probe  signals  are  seen  to  evolve  more  gradually.  The  quasi-periodic  bursts 
correspond  to  rapid  radial  transport  of  electrons  from  the  belt  region.  This  is  observed 
directly  with  a  gridded  electron  detector  and  indirectly  by  large  and  negative  decrease  in 
the  probe’s  floating  potential. 

By  using  very  high-speed  data  recorders,  we  have  identified  the  drift-resonant 
instabilities  inducing  this  rapid  electron  transport.  Figure  3  shows  a  frequency 
spectrogram  of  the  instabilities  occurring  both  during  the  microwave  heating  and  during 
the  afterglow.  During  the  heating,  the  quasi-periodic  pulsations  consist  of  relatively 
wide-band  signals  ranging  from  0.1  MHz  S/<  2  MHz.  In  contrast,  the  instabilities 
observed  during  the  afterglow  consist  of  a  multi-mode  collection  of  relatively  coherent 
rising  tones.  Closer  examination  of  the  quasi-periodic  bursts  during  heating  also  show  a 
collection  of  rising  tones.  The  difference  in  the  frequency  spectra  can  be  linked  to  the 
energy  of  the  electron  belt  as  measured  with  the  x-ray  proportional  counter.  When  the 
average  belt  energy  is  relatively  low,  (E)  -  10  keV,  (such  as  found  during  the  heating), 
the  wave  frequency  is  also  relatively  low.  During  the  afterglow,  when  the  cooler  electrons 
scattering  into  the  terrella’s  poltir  regions,  (E)  increases,  and  the  frequency  of  drift- 
resonant  instabilities  also  increase.  The  rate  of  rise  of  the  frequency  is  also  linked  to  the 
average  energy  of  the  trapped  electron  belt. 

Multiple  probes  are  used  to  determine  the  azimuthal  and  radial  structure  of  the  drift- 
resonant  waves.  When  the  probes  are  located  on  the  same  flux  surface,  the  two  probes 
indicate  that  the  waves  propagate  in  the  electron  drift  direction,  and  the  wave  spectrum 
consists  of  multiple  azimuthal  mode  numbers,  m,  as  well  as  multiple  frequencies.  When 
the  probes  are  separated  radially,  we  observe  the  radial  mode  structure  to  be  relatively 
broad  with  complex  pha.se-fronts  characteristics. 

Although  the  CTX  experiment  has  been  operating  for  less  than  one  year,  our  early 
results  already  have  important  implications  for  understanding  the  basic  process  of 
collisionless  radial  transport  in  a  dipole  magnetic  field: 

•  First,  transport  in  a  dipole  magnetic  field  seems  to  require  multiple  modes. 
This  is  a  unique  property  of  the  dipole  magnetic  geometry  resulting  from  the 
strong  radial  dependence  of  a  particle’s  azimuthal  drift  frequency,  cod  L~^, 
where  L  is  the  equatorial  radius  of  a  flux  surface.  Global  transport  can  only 
occur  with  multiple  modes  since  a  particle  does  not  remain  correlated  with  a 
single  wave  as  it  diffuses  radially. 

•  Secondly,  the  radial  transport  rate  is  fast.  In  CTX,  large  electron  bursts 
lasting  only  lO’s  of  drift  periods  are  observed  to  cause  significant  transport. 
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Figure  1.  Data  from  a  typical  pulsed  discharge  of  the  CTX  experiment 
during  the  production  of  an  energetic  electron  belt  using  microwave 
heating,  (a)  Langmuir  probe  ion  saturation  current,  (b)  hydrogen  gas 
pressure,  (c)  reflected  microwave  power,  (d)  forward  power,  (e)  probe 
fluctuations,  (0  x-ray  spectra  measured  every  50  msec. 


.5500  .5700  .5900  .6100  .B3CC  .6500 


Time  (sec) 


Figure  2.  A  fast-time  scale  detail  of  the  fluctuations  of  the  ion  saturation 
current  during  the  heating  and  during  the  "aficrglow". 
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Figure  3.  Frequency  spectrograms  of  rising-tone,  multiple-mode  drift- 
resonant  instabilities.  The  drift-frequency  of  the  energetic  electron  belt  is 
approximately  1  MHz,  and  multiple-probe  measurements  show  the  waves 
propagate  in  the  electron  drift  direction. 
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•  Finally,  drift-resonant  instabilities  lead  to  transport  causing  self-stabilization 
and  continuous,  quasi-periodic  bursting.  This  is  observation  is  important.  It 
implies  that  the  transport-inducing  waves  also  decrease  the  radial  pressure 
gradient  of  the  electron  belt — but  do  not  destroy  electron  confinement 
altogether.  For  a  dipole  magnetic  field,  the  marginally  stable  pressure  profile 
scales  like  p 

THEORY  RESEARCH  ON  SUBSTORMS 

During  substorms,  the  magnetic  configuration  of  the  magnetospheric  plasma  is 
rearranged  drastically.  The  large  electromagnetic  disturbances  and  fluxes  of  energetic 
particles  can  interfere  seriously  with  the  performance  of  spacecrafts.  As  discussed  in  our 
original  proposal  (Section  2.1),  the  ion  tearing  instability  has  been  considered  by  many  as 
a  trigger  for  substorms  over  the  last  two  decades.  We  have  treated  this  problem 
analytically,  extending  and  questioning  earlier  results  significantly. 

We  have  revisited  the  problem  of  collisionless  tearing  in  the  earth’s  magnetotail, 
and  gone  beyond  extant  theories  in  two  significant  ways.  Firstly,  we  have  included  the  By 
field  allowed  for  a  spatially  dependent  Bn  in  our  equilibrium  model.  Secondly,  we  have 
demonstrated  that  the  ion  tearing  instability,  which  has  been  suggested  by  many  as  a 
possible  trigger  for  substorms,  does  not  exist.  W'e  find,  instead,  that  if  there  is  any  form  of 
collisionless  tearing,  it  must  be  the  electron  tearing  instability.  We  have  delineated 
regimes  in  which  the  electron  tearing  mode  can  be  excited.  The  conditions  that  favor  the 
growth  of  electron  tearing  are  those  in  which  the  By  field  is  large  and  the  Bn  field  is  very 
small.  However,  if  the  Bn  field  is  significant,  the  electron  tearing  instability  is  weak,  can 
account  for  the  slow  growth  phase,  but  not  for  the  rapid  current  disruption  and 
dipolarization  phase. 

The  results  have  significant  implications  for  observations.  They  demonstrate  that 
the  By  field  must  be  included  in  a  global  model  of  magnetotail  equilibria,  and  call  for 
systerfiatic  studies  correlating  By  and  Bn  fields  with  the  occurrence  of  substorm  onset. 
Our  results  also  question  the  prevalent  wisdom  among  many  that  the  collisionless  tearing 
instability  can  provide  a  universal  trigger  for  substorms. 

For  further  details,  we  refer  to  the  attached  paper  which  is  to  be  published  in  the 
Journal  of  Geophysical  Research. 


Global  Asymptotic  Equilibria  and  Collisionless  Tearing  Stability 

of  Magnetotail  Plasmas 

XIAOGANG  WANG  AND  A.  BHATTACHARJEE 

Dtpartmtiu  ^ Appiitd  Phytia,  Cottonhia  Univtrtity,  New  York,  New  York 


Afympcodc  tail  e^oOibria  which  ate  ilowly  varying  in  the  Eanh-Son  ditechon  are  conilnioed,  including 
all  thiee  cooiponenti  of  the  magnetic  field.  Tbete  equilibria  allow  for  ipatial  dependencicj  in  and  By  . 
Some  of  these  equilibria  are  analyzed  for  ».ibiliiy  arith  reipea  to  colliiionlest  electron  tearing  moidei  using  a 
fluid  model  which  predicu,  to  within  a  nomerical  factor  of  the  growth  rates  derived  from  kinetic  theory. 
No  ion  tearing  insubility  is  found.  The  By  field  is  demonstrated  to  have  a  destabilizing  effect  on  electron 
tearing  modes.  In  the  asymptotic  equilibria  considered  here,  electron  tearing  modes  can  grow  in  the  presence 
of  By  in  those  regions  where  the  subilizing  efiect  of  electron  bounce  is  small.  Implications  for  numerical 
simulations  and  otMervations  ate  discussed. 


1.  IjflllODUCnON 

Ever  since  Ness  [1965]  reported  observational  evidence  for  a 
neutral  sheet  in  the  Earth's  magnetotail,  the  collisionless  tearing 
instability  has  claimed  much  attention  as  a  possible  mechanism 
for  magnetic  recormection  in  the  taU.  Coppi  et  at.  [1966] 
considered  a  simple  neutral  sheet  in  which  oppositely  directed 
magnetic  fields  Bj  =  (z)  face  each  other  across  the  z  =  0  line 

(in  the  x  -  z  plane)  and  demoiutrated  that  such  a  sheet  is 
unsuble  to  the  collisionless  tearing  instability  [Furth,  1962; 
Laval  et  al.,  1966].  (We  use  here  the  standard  solar 
magnetospheric  coordinates  (x,y,z),  with  the  x  axis  in  the 
Earth-Sun  direction,  the  z  axis  in  the  south-north  direction,  and 
the  y  axis  ,  which  defines  an  ignorable  direction,  is  chosen  to 
make  the  coordinate  system  right-handed.)  In  the  simple  magnedc 
geometry  considered  by  Coppi  et  al.,  the  neutral  line  z  =  0  is  the 
source  of  the  separatrix.  Far  away  from  the  separatrix.  the 
plasma  obeys  the  ideal  magnetohydrodynamic  (MHD)  equations. 
The  departure  from  ideal  MHD  behavior  occurs  in  a  narrow 
region  near  the  separatrix.  By  considering  the  energetics  of  the 
instability,  Coppi  et  al.  demonstrated  that  the  dominant 
contribution  to  the  inverse  Landau  damping  effect  comes  from 
electrons,  not  ions.  Thus  this  mode  came  to  be  known  as  the 
“electron  tearing”  instability. 

We  consider  now  the  effect  of  a  large,  constant  By  field 
superimposed  on  the  model  of  Coppi  et  al.  The  presence  of  By 
introduces  magnetic  shear  in  the  model.  The  separatrix  in  the  x  - 
z  plane  grows  out  of  the  z  s  0  line.  Drake  and  Lee  [1977] 
showed  that  collisionless  tearing  modes  are  unstable  in  this 
geometry.  (Strictly  speaking,  the  results  of  Drake  and  Lee  [1977] 
hold  for  a  low-beta  plasma  without  temperature  gradients,  as 
shown  by  Cowley  et  al.  [1986]).  Since  the  electroru  carry  the 
perturbed  parallel  current  near  the  separatrix  and  provide  the 
mechanism  for  reconnection  through  their  small  but  finite  inertia, 
this  instability  too  can  be  classified  as  electron  tearing. 

Neither  of  the  two  cases  discussed  above  are  represenudve  of 
the  Earth's  magnetotail.  Much  emphasis  has  been  placed  in  the 
literature  on  the  two-dimensional  model 


with  constant  Sq,  X,  and  [Schindler,  1974;  Galeev  and 
Zelenyi,  1976;  Lembege  and  Pellat,  1982;  Buchner  and  Zelenyi, 
1987;  Buchner  et  ai,  1991;  Pellat  et  al.,  1991;  Kuznetsova  and 
Zelenyi,  1991).  For  nonzero  values  of  B^  ,  this  configuration 
has  no  magnetic  sepaiatrix.  (For  useful  discussions  of  the  role  of 
separatrices  in  magnetic  reconnection,  we  refer  the  reader  to 
Greene  [1988]  and  Lou  and  Finn  [1990].)  The  absence  of  a 
separauix  implies  that  recormection  or  tearing  (in  the  sense  of 
affecung  a  topological  change)  cannot  really  happen  for 
significant  values  of  It  is  widely  believed  that  the  *ion 
tearing'*  instability  can  occur  in  these  circumstances,  but  the 
subject  remains  a  matter  of  lively  debate  [Schindler,  1974;  Galeev 
and  Zelenyi,  1976;  Coroniti,  1980;  Lembege  and  Pellat,  1982; 
Buchner  and  Zelenyi,  1987;  Pellat  et  al.,  1991;  Kuznetsova  and 
Zelenyi,  1991].  The  analysis  given  in  this  paper  oims  oat  to 
support  the  point  of  view  recently  expressed  by  Pella  et  al. 
[1991],  who  have  questioned  the  existence  of  the  ion  tearing 
mode.  This  point  of  view  has  significant  implications  for 
electromagnetic  particle  simulations  of  collisionless  tearing 
[Xerasavra,  1981;  Hamilton  and  Eastwood,  1982;  Swift,  1983; 
Ambrosiaru)  et  al.,  1986;  Swift  and  Allen,  1987;  Pritchett  et  al., 
1989;  Zwingmann  et  al.,  1990;  Pritchett  et  al.,  1991]  which  can 
shed  valuable  light  on  this  controversial  issue. 

Recently,  we  showed  that  the  inclusion  of  a  constant  Bj  field 
in  the  model  (1)  can  qualitatively  change  its  stability  properties 
[Wang  et  al.,  1990;  hereafter  WBL].  Our  model,  referred  to  here 
as  the  three-component  model,  breaks  with  the  tradition  of  using 
two-component  models  in  theoretical  analyses  of  collisiooless 
instabilities  in  the  magnetotail.  However,  just  as  in  the  two- 
component  model,  the  formation  of  a  magnetic  separatrix  is 
itthibited  by  the  presence  of  a  significant  B  field.  It  is 
therefore  no:  iurprisLi*  -hat  we  found  tha:  -he  electron  ’xaring 
instability  has  a  stabler  parameter  space  and  is  much  harder  to 
excite  in  the  three-component  magnetotail  than  in  a  configuration 
with  5,  =  0  .  Furthermore,  the  growth  rate  of  the  mode  is  slow, 
consistent  with  the  growth  phase,  but  not  the  expansion  phase  of 
a  substorm. 

We  now  develop  an  equili'oiium  model  which  is  more  realistic 
than  considered  heretofore  in  analytical  studies  of  collisionless 
tearing.  As  in  WBL,  we  include  the  field.  Fairfield  and 
others  have  noted  that  the  Bj  field  is  a  persistent  feature  of  the 
magnetotail  [Fairfield,  1979;  Catteil  and  Mozer,  1982;  Lui, 
1983;  McComas  et  al.,  1986;  Tsurutani  et  al.,  1984;  Sibeck  et  al., 
198S].  Voigt  and  coworkers  [Voigt  and  Hilmer,  1987;  Hau  and 
Voigt,  1992]  have  shown,  based  on  some  analytical  and 
numerical  examples,  that  the  requirements  of  global  magneiosuiic 
equilibrium  of  the  magnetotail  should  include  By  The  possible 
role  of  the  By  field  in  observations  of  substorm  dynamics  was 
pointed  out  as  early  as  1978  by  Akasofu  and  coworkers 
[Akasofu  et  al.,  1978],  From  a  detailed  examination  of  IMP  data, 
Akasofu  et  al.  stressed  the  need  for  a  three-component  analysis  of 
magnetic  fields  in  the  magnetotail  during  substorms.  Though 
there  are  a  number  of  other  events  reported  in  the  literaDrc  in 
which  an  enhancement  in  By  and  a  reduction  in  8,  is  observed 
prior  to  the  onset  of  a  substorm  [Nishida  et  al.,  1983;  Bieber  et 
al.,  1984;  Lepping,  1987;  TaJeahashi  et  al.,  1987;  Lui  et  al., 
1988;  l.ope:  et  al.,  1989],  no  systematic  studies  of  substorm 


events  with  correlated  variations  in  By  and  are  available 
yet. 

One  of  the  main  improvements  of  the  present  paper  over  WBL 
is  the  development  of  asymptotic  equilibria  with  spatial 
dependencies  in  S,  and  By ,  that  is,  5,  =  fl,  (jc.r)  and  21^  = 
By  (x.z)  .  These  equilibria,  and  the  single-particle  motions  in 
them,  are  described  in  section  2.  We  show  that  these  equilibria 
change  qualitatively  our  current  understanding  of  coUisionless 
instabilities  in  the  magnetotail  by  introducing  new  global 
features  not  captured  adequauly  by  the  model 

B  =  B.  tanh^  x  +  B  y  +  B  r  .  (2) 

with  By  and  B„  constinL  In  this  geometry,  the  global  bounce 
period  for  electrons  is  much  shorter  than  the  growth  tune  of 
the  instability  when  By-B,  .  WBL  has  been  criticized  for 
neglecting  the  subilizing  effect  of  this  bounce  [Buchner  et  al., 
1991;  Pritchett  el  al.,  1991].  This  criticism  would  be  justified 
except  for  the  fact  that  the  model  (2)  is  itself  globally  rather  crude 
and  underestimates  severely  the  bounce  period  in  a  stretched 
magnetotail.  If  we  must  include  the  effect  of  election  bounce,  then 
it  is  preferable  to  do  so  in  an  equilibrium  model  which  capcnres 
the  global  features  of  the  magnetotail  with  greater  realism  than 
equation  (2).  And  that  is  precisely  what  we  achieve  by  allowing 
for  spatially  varying  B,  in  the  new  equilibrium  model.  We  then 
show  that  there  ue  regions  where  yx^,  2  1  .  *nd  the  electron 
tearing  mode  grows,  with  a  growth  rate  y,  as  predicted  by  WBL. 
Furthermore,  when  we  coruider  the  special  case  B.  «  cottst  in 
which  yte, «  1  ,  we  find  no  ion  tearing,  contrary  to  the 
findings  of  Buchner  et  al.  [1991].  The  persistence  of  slow 
electron  tearing,  and  the  absence  of  ion  tearing,  are  recurrent 
themes  that  are  explored  in  detail  for  both  two-component  (B^  = 
0)  and  three-component  {By  ^  0)  equilibria  in  sections  4  and  3. 

Though  the  electron  tearing  instability  can  account  for  the 
growth  phase  of  a  substorm,  it  is  not  sufficiently  rapid  to 
account  for  the  cuireiU  disruption  and  diversion  that  occun  at  the 
onset  of  the  expansion  phase  \Xohahashi  et  al.,  1987;  Lui  et  al^ 
1988].  Elsewhere  [Wang  et  al..  1991],  we  have  discussed  that 
nonlinear  mode-coupling  effects  may  lead  to  a  significant 
enhaiKement  of  the  linear  growth  rate.  It  was  implicitly  assumed 
in  that  discussion  that  a  linearly  unsuble  mode  will  grow  to 
sufficiently  large  amplitude  that  it  can  couple  to  other  unstable 
modes.  However,  that  possibility  was  explored  in  the  context  of 
the  equilibrium  represented  by  (2)  which,  as  discussed  here,  has 
certain  limitations.  We  suggest  here  that  the  ciicumstaiKcs  in 
which  collisionless  tearing  can  grow  as  a  robust  iiutability.  not 
only  in  the  linear  regime  but  also  in  the  nonlinear  regime,  must 
involve  reduction  of  B,  to  extremely  small  values.  We  therefore 
point  to  the  possibilities  inherent  in  the  three-component 
asymptotic  equilibria  calculated  in  this  paper.  An  interesting 
property  of  these  equilibria  is  that  B,  may  be  reduced  to  zao  at 
near-Earth  distances  when  By  is  space-dependent.  This 
possibility  wu  first  noted  by  Hau  and  Voigt  [1992]  whose 
profiles  for  By  were  different  from  ours,  and  who  concluded  that 
for  their  class  of  profiles,  the  By  value  required  to  reduce  B,  to 
zero  is  much  larger  than  the  average  value  of  By  in  the  plasma 
sheet.  For  our  class  of  profiles,  we  find  that  the  reduction  of  B, 
to  zero  occurs  for  values  of  By  more  in  accord  with  observed 
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values.  When  5,  vanishes,  a  separatrix  can  be  formed,  and 
tearing  iiutabilities  can  grow.  Whether  the  nonlinear  evolution 
of  these  instabilities  can  actually  account  for  the  rapid  current 
disruption  and  diversion  observed  in  near-Earth  regions  remains 
an  open  question. 

The  layout  of  this  paper  is  as  follows.  In  section  2  we  obtain 
some  asymptotic  tail  equilibria,  both  with  and  without  By .  In 
section  3  we  develop  a  fluid  model,  including  a  generalized 
Ohm's  law,  which  allows  the  treatment  of  collisionless  tearing 
modes,  and  benchmark  the  predictions  of  this  model  with  known 
results  from  kinetic  theory.  In  section  4  we  use  the  energy 
integral  derived  from  our  fluid  model  to  analyze  the  stability  of 
the  two-dimensional  magnetotail  without  By  and  constant  B„ . 
We  And  that  the  ion  tearing  mode  does  not  occur,  and  the  only 
possible  instability,  under  certain  conditions,  is  an  electron 
tearing  mode.  In  section  5  it  is  shown  that  the  inclusion  of  By 
and  a  spatially  varying  B„  can  further  destabilize  the  electron 
tearing  instability.  We  conclude  in  section  6  with  a  summary  of 
OUT  results  and  their  implications  for  observadons. 

2.  Some  Asymptotic  Equiubiua 

2.1.  Two-Componeni  Equilibria 

We  begin  by  considering  some  asymptotic  two-component 
equilibria,  i.e.,  equilibria  for  which  By-Q  .  Assuming  that  the  y 
coordinate  is  ignorable,  the  magnetic  Qeld  B  may  be  written  as 

B  =  Vyxy  .  (3) 

where  V  is  a  flux  fuiKUon.  In  equilibrium,  for  a  charged  panicle 
of  type  a .  the  energy  //<,  and  the  momentum  P^y  are  conserved. 
Assuming  that  there  is  no  equilibrium  electrostatic  field,  these 
constants  are 

«a  =  2'”a^a-  W 


where  is  the  mass  of  the  charge  particle  of  type  a  ,  ,  its 

charge,  and  c  is  the  speed  of  light.  An  equilibrium  distribuuon 
funcQon  can  be  written  as 

fa  =  fa  (Ha- V  ' 


where  is  the  drift  velocity,  Tq  is  the  temperature  (in  energy 
units)  for  particles  of  type  a  {  =  e,i  for  a  hydrogen  plasma),  and 
no  is  the  average  number  density  of  electrons  and  protons.  The 
temperature  Tq  is  taken  to  be  a  constant.  From  the  reladon 

'•a  =  f‘'Va' 

and  the  requirement  of  quasi-neutrality,  n,  =  n, ,  we  obtain  the 
condition  /T,  =  -  V,y/T, .  The  y  component  of  Ampere  s 
law  gives 


(8) 


=  - —  S  ?„  J  v  /„  . 

c„  oJ  a  }■' a 


where 


(9) 


defines  a  characteristic  equilibrium  length  scale,  and  the  constant 
Bq  is  determined  by  the  relation 


=  n,  (T.  +  rj  .  (10) 


So  is  a  measure  of  the  lobe  magnetic  field.  We  scale  the  variables 
xA  -♦  X.  V/Bq  X  -»  V|/  ,  pliBoV^n)  -*  p  ,  B/Bq  -♦  B  .  and 
introduce  a  small,  positive  parameter  e  -  d/3x  «  d/3r  -  I  .  In 
this  approximation,  a  large  class  of  equilibrium  solutions  of  (8) 
is  given  by  (Birn  ef  a/.,  1975;  Birn,  1979;  Zwingmann  and 
Schindler,  1980;  Lembege  and  Pellat,  1982;  Zwingmann,  1983] 

V  =  -log  (cosh  [r/ez)] //(cc) }  ,  (11) 


where  flex)  U  a  slowly  varying  function  of  z  .  It  follows  that 


B 

I 


=  yiez)  tanh  [^£z)) 


(12) 


9m/  If'itx)  1 

^  («)tanh[g^ez)]|,  (13) 

where  prime  denotes  differentiation  with  respect  to  the  argument. 
Following  Lembege  and  Peliai  [1982],  we  take  /fez)  =  exp  (ez) 
(z  <  0) .  As  a  first  approximation,  almost  all  analytical  studies 
replace  (12)  and  (13)  with 


B  3  tanh  z  , 

X 

(14) 

B,=  c. 

(15) 

However,  this  is  valid  only  in  the  region  Izis  1  ,  Izl— >  0  .  It 
cannot  be  assumed  that  this  approximation  holds  for  large  z 
close  to  the  z  =  0  line,  for  the  equilibrium  pressure  balance 
condition  implies  that  dpBx  =  -  e  which,  in  turn,  gives  p  =  pt)- 
ez ,  with  po  constant.  This  means  that  the  pressure  increases 
with  z  ,  which  is  uitrealistic  for  the  distant  tail.  A  better 
approximation  for  B„  than  (15)  is 

B^  =  e  (1  -z  tanh  z) . 

(16) 

However,  (16)  does  not  represent  realistically  the 
of  B,  (zj)  which  should  tend  to  zero  as  Izl— » 
to  improve  on  (16)  by  taking 

z  dependence 
.  It  is  possible 

fitx)  =  exp  (ez/(l  -  £z)l  (z  <  0) . 

(17) 

We  recall  that  the  distance  x  is  measured  in  units  of  \~l 
For  speciftcicy,  we  take  E  a  0  1  ,  and  define  Lq  b  e~'  =  10  .  The 
scale  Lq  (-lO^c)  enables  us  to  define  three  separate  regions: 

1.  Near-Earth  region.  Ixl  «Lo 

In  this  region,  using  (17)  in  (12)  and  (13),  we  get 

-  tanhz,  (18) 

=  E  (1  -  z  tanh  z) .  (19) 

2.  Middle  region.  Ul  - 

In  this  regioa  since  EX  -  1 ,  we  get 

=  e  unh  (e  z)  ,  (20) 

fl,  =7  [l-e’^^ztenh  (e'"z)]  .  (21) 

n  4 

If  we  define  B'  a  e"*^B  ,  z*  a  and  E'  a  e/4  ,  then 
(20)  and  (21)  became 

a  ^'  =  tanh  s'.  (22) 

B,  '  =  s'  (I  -  s'  tanh  z') .  (23) 

which  is  the  same  as  (18)  and  (19).  In  other  words,  the  middle 
region  has  the  same  structure  as  the  near-Earth  region,  except  that 
the  magnetic  field  in  it  is  somewhat  weaker  and  the  current  sheet 
is  somewhat  wider. 

3.  Distant-tail  region.  lxl»Lo 

In  the  distant-tail  region,  taking  the  limit  x/Lq  we  get 

B^  =  «  '  tanh  («  *  z) ,  (24) 

B,  =  0.  (25) 

Equations  (24)  and  (25)  describe  essentially  the  configuration 
first  considered  by  Coppi  et  al.  [1966],  who  found  electron 
tearing  modes  to  be  unstable.  We  caution  that  though  the 
qualitative  features  of  the  magnetic  field  as  described  by  (24) 
and  (25)  are  reasonable,  the  far  distant-tail  region  is  outside  the 
strict  domain  of  validity  of  the  asymptotic  solution  (13).  Here 
we  do  not  pursue  this  maoer  further,  for  the  main  focus  of  this 
paper  is  on  the  investigation  of  coUisionleat  tearing  modes  in  the 
near-Earth  and  middle  regions.  We  also  note  that  the  near-Eanh 
configuration  described  above  is  not  totally  realistic  for  txl  S 
5Rg  because  at  these  distances,  the  Earth's  dipole  field,  not 
included  in  the  model,  plays  a  dominant  role.  Matching  to  the 
dipole  field  can  be  carried  out.  in  principle,  but  is  not  germane  to 
our  considerations  here.  Figure  1  shows  a  plot  of  our  asymptotic 
two-dimensional  model  with  the  origin  set  (arbitrarily)  at  z  -  - 
5  Rg  z  =  0. 

We  now  discuss  some  features  of  the  single-particle  orbits. 
For  a  particle  of  mass  mg  aitd  charge  gyrating  in  magnetic 
field  B  ,  the  Larmor  frequency  is  (Oeg  »  q^BIm^c  .  and  the 
typical  Larmor  radius  is  p^g  -  ^lal^ca  •  "'here  v,ga 
{2T is  the  thermal  velocity.  Using  typical  uil 
parameters  (see.  for  instance,  Lui  [1987]),  we  get  pj,  -  (0.5  - 


In  both  regions  1  trd  2.  pe,  2  1 ,  p„  «  1  (scaled  by  X). 
Hence  electrons  rray  be  treated  in  the  guiding-center 
approximation,  but  the  ions  are  essentially  urunagnetized. 
Because  of  the  r  dependence  of  (x.z) .  the  field  lines  are 
more  stretched  as  z  increases  than  in  the  case  of  constant  B„ 
(see  Figure  1.).  In  the  Appendix,  we  show  that  this  has  the 
consequence  that  the  avenge  bounce  period  of  electrons  in  both 
regions  2  and  3  can  be  substantially  larger  than  the  bounce 
period  with  constant  B^. 

It  is  interesting  to  note  that  the  large  separation  in  the 
magnitude  of  the  Larmor  frequency  oOg^  and  the  average  bounce 
frequency  cOs*  *  in  our  asymptotic  equilibria  (in  those 

regions  where  fr,  is  weak)  diminishes  the  possibility  of  low- 
order  resonances  between  the  Larmor  and  bounze  frequencies  for 
most  particles,  and  the  chaos  that  can  result  from  such  resonances 
[Chen  and  Palmadesso,  1986;  Bdchner  and  Zelenyi,  1987], 
Hence,  for  equilibria  with  spatially  varying  B, ,  we  will  not 
concern  ourselves  here  arith  intrinsic  stochastic  diffusion  as  a 
possible  mechanism  for  the  restoration  of  ion  tearing.  As  to 
whether  stochastic  diffusion  can  the  ion  tearing  if 

is  constant  has  been  the  subject  of  debate  recently  [Peltat  et  al^ 
1991;  KuXTieisova  and  Zelenyi,  1991]  and  is  an  issue  we  shall 
address  in  section  4. 

2.2.  Three-Component  Equilibria 

We  now  consider  equilibria  which  are  symmetric  in  y  but 
with  aU  three  components  of  the  magnetic  field  nonzero.  The 
magnetic  field  B  is  represented  as 

B  =  B^y  +  Vy  X  y  .  (26) 

The  condition  of  magnetostatic  equilibrium  gives  the  Grad- 
Shafranov  equation  [Voigt  and  Hilmer,  1987;  Paranicae  and 
Bhanacharjee,  1989;  Hau  and  Voigt,  1992] 

<27) 

where  psp(y)  and  By  ■fl,(v)  are  two  free  functions,  ffwe 
set  By  *0  and  take  p  »  exp  (2v)  .  we  recover  (8)  (in 
dimensionless  variables).  We  note  that  (8)  also  holds  for  By  <= 
const,  in  which  case  a  class  of  asymptotic  solutions  can  be  again 
consLucted  using  (11).  The  B,  and  B,  components  for  tlus 
class  of  solutions  has  already  been  given  in  section  2.1. 

We  now  consider  the  effect  of  a  nonzero  dBy  Vdy  on  the 
solution  (11).  Since  y  is  small  and  negative  near  zsO.we 
make  the  expansion 

.  (2g) 

when  ho  and  b^  are  positive  constants.  If  we  order  ho  ~  ~ 

ly  I -C ,  then  near  z  s  0  (11)  can  be  modified  as 

y  =  -  log  (cosh  [z/i:cx)]/y^cx)]  +  b^x^  11,  (29) 

with  /(ex)  specified  by  (17).  This  yields  B,  =  e-»-h|X 

which  implies  that  B,  vanishes  on  the  z  ^  0  plane  at  x  =  - 
hj/e  .  Hence  an  X-point  appears  on  the  z  »  0  line  (in  the  x  -  z 
plane)  at  near-Earth  distances  (5  Bj  <lx  I  510  Rg).  That  the 
spatial  dependence  of  B^  in  a  magnetostatic  equilibrium  can 
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lead  to  the  reduction  of  to  zero  on  the  z  »  0  plane  has  been 
recognized  by  Hau  and  Voigt  {1992].  However,  for  their  class 
of  equilibrium  profiles,  they  find  that  the  Bj  required  to  reduce 
B,  to  zero  is  much  larger  than  the  observed  average  in  the 
plasma  sheeL  For  our  class  of  pressure  and  By  profiles,  this 
limiution  is  overcome  because  the  average  values  of  both  B. 
aitd  By  are  of  the  same  order. 

The  panicle  orbits  and  drift  modons  in  the  presence  of  By  has 
been  considered  by  WBL  and  will  not  be  repeated  here.  Simple 
consideradons  of  particle  orbits  near  z  »  0  ,  where  (2)  holds, 
shows  that  the  main  stabilizing  effect  of  B,  is  to  remove 
electrons  from  the  z  =  0  plazie  where  the  separaiiix  tends  lo  form 
in  the  absence  of  B„ .  The  addidon  of  the  By  field  provides  the 
electrotu  with  another  guided  chaimel  for  modon  near  the  z  >  0 
plane.  Then,  as  noted  by  WBL.  the  growth  of  the  electzon 
tearing  instability  depends  on  the  competing  effects  of  By  tai 
B,.  As  far  as  the  effect  of  electron  bounce  is  coocemed.  we  show 
in  the  Appendix  that  the  bounce  period  in  the  stretched 
tnagnetotail  increases  due  to  the  z  dependence  of  B,  ;  hence  dte 
condidon  2  1  is  satisfied  in  the  middle  region.  This,  in 
rum.  implies  that  the  stabilizing  effect  of  the  electron  bounce  is 
weakened,  and  that  the  electron  tearing  mode  can  grow  in  the 
linear  regime  not  only  for  By  »  B,  but  also  for  By  -  B,  near  the 
z  =  0  plane.  However,  we  repeat  for  emphasis  that  in  the  latter 
case,  the  growth  rate  is  small,  and  the  mode  is  likely  to  saturate 
nonlinearly  at  a  low  amplitude. 

In  this  section  we  have  made  use  of  dimensionless  variables  in 
order  to  keep  the  noudon  simple.  In  the  remainder  of  the  paper, 
we  shall  return  to  using  the  primitive  physical  variables. 

3.  Energy  Integral 


In  WBL  the  subilicy  of  the  plasma  sheet  was  analyzed  by 
asymptouc  matching  of  the  normal  mode  equations  in  the  inner 
region,  where  finite  particle  inerda  provides  the  recoonecdon 
mechanism,  with  the  equadons  in  the  outer  regioa  governed  by 
ideal  MHO.  The  technical  details  of  such  an  approach  are 
somewhat  different  from  those  involved  in  the  Lyapunov 
funcdonal  method  (developed  by  Laval  el  ai.  [1966])  which  relies 
on  the  existence  of  an  energy  integral.  We  review,  at  first,  the 
subility  criteria  that  follow  from  the  energy  integral,  from  the 
linearized  Maxwell's  equadons,  it  follows  that 


f 


=  -  ®  da  —  E  X  B 
471  '  * 


(30) 


where  all  perturbed  quandues  are  designated  by  the  subscript  I. 
We  assume  that  the  boundary  condiuons  on  the  surface  bounding 
the  plasma  volume  ensure  that  the  surface  term  vanishes.  Then,  it 
can  be  shown  [Laval  et  al.,  1966]  that 


(31) 
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where 

df 

/la  -  Ao  '  ij  V,  . 
From  (30)  and  (31),  it  follows  that 

IjA.O. 


where 


(32) 


(33) 


The  energy  integral  S^e  is  a  quadratic  form.  If  S^e  is 
posidve-definite  for  all  nontrivial  permissible  perturbations, 
then  the  system  is  stable  [XrusJcal  and  Oberman,  \95i;  Laval  et 
al.,  1966].  In  other  words,  a  sufficient  condition  for  stability  is 

5^e  2  0  .  (35) 

Furthermore,  for  a  Maxwellian  distributkm.  since,  d/a/dffa- 
-/o/Ta ,  we  get 


5  ^ 


2  0.  (36) 


Equadon  (36)  implies  that  for  a  Maxwellian  distribution  a 
sufficient  condidon  for  stability  is 


3  dy 


2  0. 


(37) 


Since  coUisionless  tearing  modes  have  low  frequency,  the 
electrical  energy  E\^l%n  is  much  smaller  than  the  magnedc 
energy  ,  and  can  be  neglected.  The  sufficient  condition 

(37)  can  be  rewritten  as 


6^W  I 


/■  S-i??/ 


fco. 


(38) 


One  of  the  difficulties  presented  by  the  energy  integraJ  5^  is 
that  the  physical  interpretadon  of  some  terms  is  not  transparent. 
We  take  for  instance,  the  last  tetm  on  the  right-hand-side  of  (34). 
Lembege  and  Pellat  (1982)  showed,  by  using  a  Schwartz 
inequality,  that  the  term  has  a  lower  bound  which  can  be 
attributed  to  the  compressibility  of  the  electron  fluid.  This 
interpretadon  has  been  invoked  repeatedly  in  the  literature,  but  it 
is  worth  noting  that  it  was  meant  to  hold  for  (he  lower  bound,  and 
not  the  term  itself.  In  fact,  it  is  clear  that  there  is  more  to  the 
original  term  than  electron  compressibility.  The  funcdon  S^e 
represents  the  second  vahadon  of  the  total  energy  e  which  is  a 
sum  of  the  electromagnedc  field  energy  and  the  kinedc  energy  of 


the  fluid.  Since  the  perturbed  kinetic  energy  of  the  fluid  is 
always  positive  definite,  it  must  be  associated  with  a  manifestly 
positive  definite  term  in  5^  .  Thus  it  would  seem  that  the  last 
term  in  (34)  should  involve  the  perturbed  kinetic  energy  of  the 
fluid  plasma,  but  this  is  not  obvious.  (The  physical 
interpretation  of  the  kinedc  terms  is  much  more  transparent  in  the 
energy  principle  of  Kruskal  and  Oberman  [1938]  for  a  guiding- 
center  plasma,  but  this  is  not  the  model  underlying  the  functional 
(34)  which  has  been  derived  from  the  fun  Vlasov  equadon.) 

In  view  of  the  interpreudonal  dtfflculdes  of  a  fully  kinedc 
treatment,  we  propose  a  different  approach  for  the  study  of 
collisionless  tearing  modes  in  the  magnetotail.  This  approach 
uses  fluid  equations,  and  the  mechanism  for  reconnection  is 
provided  by  finite  particle  inerda  in  the  generalized  Ohm's  law. 
That  this  is  a  reasonable  model  for  coUisionleas  tearing  inodes 
was  suggested  by  Furth  [1962,  1964].  The  point  of  view  we 
adopt  here  Is  that  the  energedcs  of  collisionlesa  tearing  modes  is 
describable  within  a  fluid  model  by  using  a  generalized  Ohm's 
law.  The  wave-panicle  resonance  condition  (which  is  a  Idiietic 
effect)  is  included  in  such  a  nwdel  by  simply  equating  the  growth 
rate  y  to  kv^  •  where  k  is  the  wave  number  and  v,^  the  elecmin 
thermal  speed.  (For  sheared  configuradotts,  k  is  replaced  by  k  | , 
the  component  of  k  parallel  to  the  magnetic  field.)  This  ad  hoc 
resonance  condition  misses  the  detailed  structure  of  the  pardcle 
distribudon  funcuons.  However,  we  shall  demonstrate  that  the 
growth  rate  calculated  from  the  model  equations  agrees,  except 
for  an  overall  muluplicative  factor  of  with  the  results  of  a 
fully  kinedc  treatment.  One  of  the  advantages  of  the  fluid  model 
is  that  the  energedcs  of  the  instability  is  much  easier  to  inteqiret 
physically.  This  will  enable  us  to  formulate  a  stability  condition 
which  is  both  necessaiy  and  sufficient. 

The  linearized  fluid  equations  are 


3u  JxB,  J,  xB 

dr  c  c 


(39) 


3  P, 

+ V-(pu) 


0, 


(40) 


u  xB 


3J, 

3r 


I  dB, 


(41) 


(42) 


VxB,-^J,.  (43) 

p,  -  11,(7, +  (44) 

where  p|  is  the  perturbed  mass  density,  u  is  the  perturbed  fluid 
velocity,  pi  is  the  perturbed  pressure,  ri)  is  the  perturbed 
number  density  (/ij  s  ^  by  quasi-neutrality),  and  the 
unsubscripted  variables  represent  equilibrium  quanddes.  We 
note  that  the  generalized  Ohm’s  law  (41)  has  a  term  propordonal 
to  the  electron  mass,  but  none  propordonal  to  the  ion  mass.  This 


can  be  readily  seen  by  considering  the  more  general  form  (see,  for 
instance,  Krall  and  Trivelpiece  [1986]) 


which  reduces  in  the  limit  mjmi«  1  to  the  form  (41).  Apart 
from  terms  involving  finite  particle  inertia,  the  generalized  Olm's 
law  contaiiu  other  terms  such  as  the  Hall  term,  the  electron 
pressure  gradient  teim,  as  well  u  terms  involving  the  anisotropic 
stress  tensor.  It  can  be  shown  that  the  first  two  do  not 
qualitatively  change  our  results  for  the  class  of  equilibria 
considered  here.  However,  anisotropies  in  the  stress  tensor, 
which  are  an  additional  source  of  free  energy,  may  alter  some  of 
our  conclusions.  We  do  not  consider  pressure  anisotropies  hcR 
because  it  is  questionable  whether  an  instability  that  is  primarily 
driven  by  such  a  source  of  free  energy  should  be  classified  as  a 
tearing  mode. 

We  now  use  (39)  -  (44)  to  calculate  the  different  terms  in  (30) 
(neglecting,  of  course,  the  term  which  is  much  smaller 

titan  the  term  5i*/8Jt).  We  get 


Ji-Ei 


J,  xB  u 

c 


(43) 


The  first  term  on  the  right-hand-side  of  (45)  can  be  calculated 
from  the  momentum  equation  (39)  which  gives 


J,  xB  u 
c 


t-  V  (p|U)  -pjV  ut-u - 

(46) 

Writing 

Bj  =  VVj  X  y  -t-  y  . 

(47) 

we  get 

V, 

-  - ^  V  •  u  -H  — u  •  V/ . 

c  c 


(48) 


Since  3vi/3t  -t-  u  •  Vy  =  0  ,  and/  >  J(yf) ,  we  obtain 


u- V/ 


”  4jt  3t  v'(z.ex)  ’ 


(49) 


where  the  prime  in  \)r'(z,ex)  indicates  differentiation  with 
respect  to  z  .  Defining  F  ■  y'(z,ex)  and  substituting  (46)  - 
(49)  in  (45),  we  get 


1  2  ,  ■'t  f  1 

p 


(50) 


2 

where  s  Im,  is  the  plasma  frequency,  and 

Pi  = /»l -P'(V)  Vi  »  P,  — /¥i/c.  (51) 

Equation  (30)  can  then  be  cast  in  the  fomu 

=0.  (52) 

where 

5^e  =  5 +  5 +  S^jf  +  5^2  *  5^6r  +  S^AT  +  5^Q .  (53) 
Here 

is  the  free  energy  of  the  magnetic  Geld; 

1^  SV  -  -  J  rfxpjV-u,  (55) 

can  be  attributed  to  plasma  compression,  and 

5^2  s  Jdx(2n/(oJ)jJ.  (56) 

is  the  "dissipation"  due  to  finite  electron  inertia.  Note  that  the 
fluid  kinedc  energy  S^Af  ■  jdx  (\l2)ptP-,  and  the  dissipation 
522  are  positive  definite  quantides.  If  the  boundary  in  (52)  is 
chosen  such  that  the  surface  term  vanishes,  then  a  sufGcient 
condition  for  stability  is 

5‘C/  =  2  0 .  (57) 

In  order  to  develop  confidence  in  the  fluid  model,  we  now 
benchmark  it  with  standard  results  from  kinedc  theory.  For 
simplicity,  we  consider  an  incompressible  fluid  for  which 
B  0  .  Then,  there  are  three  terms  in  the  energy  integral  (33)  among 
which  S^K  is  the  kinetic  energy  and  always  posidve  definite.  It 
is  clear  that  a  tearing  instability  can  occur  if  and  only  if  there  is 
magnetic  free  energy  available,  i.e.,  S!^f  <  0  ,  and  there  exists 
simultaneously,  a  mechanism  for  dissipation  causing  5^ 
(which  is  always  posidve  definite). 

Let  us  first  consider  the  equilibrium  configuration  (2)  with  5, 
s  0  .  For  this  case  the  growth  rate  of  the  collisionless  tearing 
insubility  is  known  analydcally.  The  maihemaucal  problem 
divides  itself  neatly  into  two  regions:  the  outer  region,  away 
from  z  s  0  ,  where  the  plasma  is  in  quasisiadc  equilibrium,  and 
the  inner  region  near  z  «  0  where  inertial  and  dissipative  effects 
are  important.  We  rewrite  (54)  in  the  form 

where  ipj  »  vi/(z,  t)  cosJtx  ,  F  a  tanh  zA  and  prime  denotes 
derivadve  with  respect  to  z  .  In  the  outer  region,  neglecting 
inerda  and  dissipadoa  we  get 

=  0. 


(59) 


Furth  [1963]  observed  that  if  the  first  term  in  the  integrand  of 
(58)  is  integrated  by  parts  and  (59)  is  used,  reduces  to 


r*0  + 
r=0- 


“  -  Tf"  ^  ^^(0) .  (60) 

1671  ^ 

where  jdx  »  u^/<fr  ,a^  *  jdxdy  and  the  parameter  Aq' is 
defined  by  the  relation. 


Ao’  » 


v(0) 


(61) 


( In  obtaining  (60) ,  we  have  used  the  boundary  conditions 
\j/(-  m)  3  y(-f  <■)  s  0  .)  For  the  equilibrium  (2)  with  =  0  • 
get  (see,  for  instance,  WBL) 

(62) 

w  a.^  * 


On  the  interface  between  the  inner  and  outer  regions,  the 
Poynting  flux  is 


da  • 


X  B 


1 


16tt 


^  h 


L  *2 

■at  * 


(63) 


Thus  the  magnetic  free  energy  in  the  outer  region  is  spent  as 
dissipation  and  kinetic  energy  in  the  iiuier  region. 

To  determine  the  dissipation  in  the  iimer  region,  we  use  the 
generalized  Ohm's  law, 


1  ay  ^Ll 

c  dt  ~  ngi  dt  ' 


where  Jiy  *  Jy  cos  kx.  Equation  (64)  gives 


J 


J 


2 

0 


V. 


(64) 


(65) 


where  Jbo  m  co^/c  .  The  dissipation  caused  by  electron  inertia  in 
the  inner  region  is 

-  r=+4, 

J  *;  r  -I 

TT.  J 

r 


20 


*0  -  2 

=  ^  V(0)  .  (66) 

where  d,  is  the  width  of  the  tearing  layer. 

Requiring  that  the  system  energy  be  equal  to  its  equilibrium 
value  [Kruskal  and  Oberman,  19S8],  we  get 

=  sVy.  +  +  8^2  =  0 .  (67) 

For  the  electron  tearing  instability,  the  fluid  kinetic  energy  f?’K , 
which  is  dominantly  due  to  ions,  is  much  smaller  than  the  electron 
dissipation  5^ .  Hence,  from  the  relation 

+5*2  =0,  (68) 

we  obtain  the  tearing  layer  width  [DraJte  and  Lee,  1977] 

4 

d  =  d(k)  =  (69) 

2*0 

The  growth  rate  of  the  instability  can  be  determined  by  using 
the  wave-particle  resonance  condition.  In  the  case  Bj  =  Q  this 
condition  is 


y  =  kv  , 

where  k  is  determined  from  (62)  by  the  requirement 
(S^Wy  <  0).  From  (62)  and  (70)  we  obtain 


—  (i-*V). 


1/2  / 

I 


T, 

T 


l-i-=r-  (1-xV). 


(70) 

Ao'>0 

(71) 


(72) 


In  writing  (72)  we  have  made  use  of  the  equilibrium  relation  (10). 
The  growth  rate  (72)  is  larger  by  the  factor  than  the  result 
obtained  from  kinetic  theory  by  Laval  et  al.  [1966]. 

In  the  case  By*Q  the  resonance  condition  (70)  is  modifted  to 


7  = 


(73) 


where  Jkii  =  kd,[k  .  Equations  (62)  and  (73)  then  yield 


T  = 


2t*X 


(74) 


which  is  also  larger  by  the  factor  than  the  result  of  Drake 
and  Lee  [1977]. 

Equations  (71)  and  (74)  demonstrate  that  the  fluid  model  is 
reliable  as  a  predictor  of  the  parametric  dependencies  and  the 
order  of  magnitude  of  the  growth  rate  of  collisionless  tearing 
modes  when  B,<=0  ,  both  with  and  without  By .  The  numerical 
factor  missed  by  the  model  has  to  do  with  the  precise  details  of 
the  electron  distribution  function. 

In  the  last  pvagraph.  we  have  used  the  term  "collisionless 
tearing"  instead  of  "electron  tearing."  though  the  latter  name  is 
commonly  used.  It  has  been  shown  elsewhere  (for  the  case  By  * 


0)  that  Che  full  ion  dyntmica  gives  small  logarithmic  conections 
to  (74)  (WBL.  Appendix  A].  In  other  words,  the  result  (74) 
already  includes  the  ion  response,  and  there  is  no  ion  tearing 
branch  of  the  dispersion  equation.  That  being  so.  it  is  redundant 
to  call  this  tearing  instability  electron  tearing  because  in  this  case 
there  is  no  other  form  of  tearing. 

4.  Effect  of  Constanta, 


So  far.  we  have  benchmarked  the  fluid  model  by  reproducing 
known  results  from  kinetic  theory  for  equilibria  with  A,  =  0  . 
We  presume  that  a  fluid  model  whidi  hu  been  demonstrated  to  be 
reliable  for  A,=  0  will  work  for  equilibria  with  A,*  0. 
representative  of  the  Earth's  magnetotail.  As  discussed  in  WBL, 
what  will  change  in  the  presence  of  a  nonzero  A,  are  the  particle 
orbits.  We  consider,  at  first,  the  case  A,  ^  cotutant  (equation 
(15))  which  has  engendered  coruiderable  controversy  in  the 
literature. 

In  the  presence  of  A,  the  perturbed  current  has  a  component 
7(1  in  the  r  direction  which  generates  a  perturbed  nugnetic 
field  component  Biy.  In  WBL,  this  component  hu  been 
calculated  by  writing  Bj  s  V  x  Ai ,  using  the  Coulomb  gauge  V 
*  A 1  =  0  and  the  approximation  d/3x  «  d/3z  for  perturbed 
quantities.  Herewedenote  yisAiy,  Xi 


Ai  a  Vl  + 


(75) 


Substituting  (75)  in  (54),  we  have 


SV  =  J  dx 

I  8rt  ■' 


.2 

V,  + 


f"  2  .2 


(76) 


We  take  xj  =  Xi(z.t)  cos  tx  .  In  the  outer  region,  we  have 
[WBL,  Appendix  B] 


x"-k^X=0.  (77) 

which  hu  a  solution  of  the  form  .  This  solution  hu  a  jump 
continuity  in  its  logarithmic  derivative,  specified  by 


^2“  - 1 - 

X(0) 


-2k. 


(78) 


In^the  Coulomb  gauge,  we  can  rww  write  Aj^*  =  (X  -k 
X)  .  Then,  the  contribution  of  the  lut  term  in  the  integrand  of 
(76)  is  -  {Ojy  /16r)  ^2  X^(.0)  .  This  hu  a  clear  physical 
iruerpretation;  A,  reducu  the  magnetic  free  energy  available  to 
the  tearing  instability  and  hence  hu  a  stabilizing  effect. 

The  computation  of  the  first  three  letiru  in  the  integrand  of  (76) 
is  somewhat  more  involved.  One  of  the  complications  introduced 
by  A,  is  that  it  introducu  a  phase  shift  that,  in  effect,  couplu 
the  cosine  solution  in  x,  i.e.,  (z.r)  cos  kx  with  the  sine 

solutiotu  in  x,  i.e.,  y,  (rx)  sin  )bc .  In  the  outer  region  equation. 


y'l^v'yj  =y" 


3y, 

dx 


subsiiiuung 


(79) 
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coj«  im  i= 


(80) 


we  obtain 


A"  v_.  a  A'  (w_  -  i't/^  " 


A”'  a  A'(v  -i'v)  -  t'vp.  (82) 

where  A'  ■  X~‘  tanhz  /  X  .  Fotlowing  the  aethod  outtinad  is 
Appendix  3  of  WBL  we  get 

=  v;®’  ±  cXAA'^^®’  +  0(0 .  (83) 


=  55;°’  =  eXAA'vI*”  1-0(0.  (8^) 

where  the  upper  sign  in  (83)  and  (8d)  ccneuxinda  va  z  ^  0 ,  ±e 

lower  sizn  to  z  <  0  ,  and  b  a  d/(A  d^'  \?) .  The  leadine  oicer 
,  .  -.0)  .  '  '  * 
somaons  and  ooey 

A  '  Wjj  =  A  .  (82) 

We  now  use  (81)-  (82)  a  calculate  ie  fiia  three  amis  in 
.  These  are 


5"W  3  ”  I  ax  r  'v.  •■ 

Zr,  J  L  • 


M _ : 


»?]■ 


.2-*  jr.i 

* -  ax  V,  +•<  p  1- — p, 

16k-  «  A' 

,2.1  A-.zl 
-  v,  -  i  • 


(86) 


Integrating  the  right-hand  side  of  (86)  by  parts  and  using  (81) 
and  (82),  we  get 


a  ^  .  r  rr  5'“’ 
51V  r--^  d,'‘5*(0)  lldz^ 
»  Idrr  ^  a  U  V 

r.  «”  /  V 


W'”’ 

K  • 


(87) 


where  p  (0)  m  \5j*(0)  ”  9,^(0) .  Since  p,(®'  and  V,  ®  ire 
both  solutions  ot  (82)  and  obey  the  same  bonndary  coodidonx. 
the  last  two  terms  of  (37)  cancel  each  otfas  exactly.  DeSning 


d,  a 


and  using  (87),  we  get 


A,  7^0)  ■»-d.£*(0) 
7  (0)  l‘(0) 


(88) 


5'-''^  r  -■—  d'  (p^Ol^UO))  .  (89) 


We  now  evaluate  ^  ,  which  is  due  to  plasma 

compressibility.  If  is  large  and  constant,  then  the  inequality 
«  1  is  satisfied.  Under  these  conditions  the  electrons  have 
a  stabilizing  compressional  contribution  due  to  the  bouncing 
motion  between  mirror  points  along  a  field  line.  However,  the 
effect  of  ion  compressibility,  which  comes  into  play  because  of 
the  quasi-neutrality  constraint  =  /i,|  ,  is  larger  than  the  effect 
of  electron  compressibility.  The  magnetic  free  energy  5^^  is 
generally  not  large  enough  to  provide  energy  for  compressing  the 
ions,  and  the  mode  is  subilized  unless  the  wavelength  is  very 
large.  In  order  to  demonstrate  this,  we  calculate  from  the 

fluid  equations.  From  the  linearized  continuity  equation,  we  get 


1  d 

V-u  =  --  -  u  •  vw  —  Inn  . 

n  dt  dy 

Using  the  equilibrium  relation  (7),  and  defining 
_  dn 

't  -  — 

we  rewrite  (90)  in  the  form 


n.-—  Vj. 


(90) 


(91) 


1  3n  2  1^1  _  I 


(92) 


Averaging  (92)  over  an  electron  bounce  period,  we  have 


<^-u>  = <*‘U>)  .  (93) 

where  <  >  indicates  an  average  over  the  rapid  bounce  motion  of 
the  electrons.  Since  the  bounce  modon  involves  dominantly  the 
outer  region,  we  neglect  the  electron  inertia  term  in  averaging 
Ohm's  law  (41).  We  then  obtain 


<E  > 


c  ^  dt  ^  - 


(94) 


Using  (44).  (51),  (93),  and  (94),  we  get 

<?,  v.u>  =  <1^  vf> 


16ii  B 


<3  2 


(95) 


Hereafter,  to  simplify  notadon.  we  shall  drop  the  averaging  sign. 
Equadon  (95)  then  gives 


dx 

16n 


2 


-2 
V  • 


(96) 


Equadon  (96)  can  be  rewritten  as 


1 6rt  2fl  ^ 

fl 


to  V^O) . 


(97) 


where 


(98) 


Zq  =  2  j  dz  exp(-2kz)  [  1  +  tanh  zA]  . 

Equation  (97)  agreu  with  (91)  of  GaUev  [1984]  (except  for  the 
area  term  a^y  which  haa  been  taken  to  be  unity  by  Galeev 
without  loss  of  generality).  The  constant  zq  is  estimated  by 
Galeev  from  physical  arguments;  here,  we  evaluate  (98) 
asymptotically  to  obtain 

*0  =  lA  .  (99) 

Note  that  an  upper  bound  for  is  obtained  by  setting  % 

(0)  equal  to  zero.  Hence,  from  (57),  a  sufficient  condidoa  for 
stability  is 

-  Ag'  +  kjz^  i.  0  ,  (100) 

which  for  £  «  1  ,  reduces  to 

kX/Z  >2.  (101) 

Equation  (101)  is  close  to  the  sufficient  condition  for  stability 
kX/z  2  4/jt  ,  derived  by  Lembtge  and  Pellat  [1982].  The 
inequality  (101)  implies  that  all  wavelengths  smaller  than  itXjz 
are  suble.  For  example,  if  we  take  X.  -  1  ^£,  £  -  0.1  ,  we  find 
that  wavelengths  smaller  than  30  R[  uo  stable.  Of  course,  this 
does  not  necessarily  mean  that  wavelengths  larger  than  irX/z  are 
unsubte  because  violation  of  (101)  does  not  imply  insubOity.  If 
one  proceeds  with  the  hypothesis  that  instability  is  possible  for 
kX/z  <  2  ,  it  can  be  shown,  following  Lembtge  and  Pellat  (1982], 
that  a  long-wavelength  ion  tearing  mode  is  impossible.  We  refer 
the  reader  to  the  work  of  Lembege  and  Pellat  [1982]  for  further 
details. 

Anempts  have  been  made  to  restore  the  ion  tearing  insubOity 
by  invoking  pitch-angle  diffusion  [Coroniti,  1980;  Galeev,  1984 
and  references  therein]  or  intrinsic  chaotic  diffvuion  [Buchner  et 
al.,  1987].  We  now  demonstrate  that  even  in  the  presence  of 
these  effects,  the  most  that  we  can  get  is  some  form  of  weak 
electron  tearing  and  that  there  is  no  ion  tearing.  At  first  we  note 
that  occasionally,  a  source  of  some  confusion  in  the  literature  has 
been  the  misleading  premise  that  it  is  electron  compressibility 
that  stabOizes  tearing  in  the  presence  of  a  constant  .  From 
this  premise  follows  the  argument  that  if  the  electrons  are 
removed  by  pitch  angle  scattering  or  intrinsic  stochastic 
diffusion,  then  it  is  possible  to  neglect  the  electrons  while  the 
ions  tear  field-lines.  (3ur  fluid  model  clearly  indicates  that 
electron  compressibility  is  less  of  a  factor  than  ion 
compressibility  for  conditions  typical  of  the  magnetotail. 
Inspection  of  (95)  shows  that  both  ioru  and  electrons  contribute 
to  ,  but  the  electron  contribution  to  may  be 

apportioned  as  [T, /(T.-f  T,]  ^  ,  whereas  the  ion 

contribution  is  T,\  .  Since  7,  ■  57,  is  typical  in 

the  magnetotail,  this  apportionment  indicates  that  the  dominant 
contribution  to  fluid  compressibility  comes  from  ions. 

In  order  to  pinpoint  the  differences  between  our  results  and 
others  in  the  literature,  we  refer  the  reader  to  the  review  by 
Galeev  [1984],  Galeev's  equation  (91)  gives  the  energy  spent 
for  pluma  compression,  in  agreement  with  our  . 

Subsequently,  in  the  presence  of  pitch  angle  scattering,  Galeev 


■ctributes  a  compressional  term  similar  to  entirely  to 

electrons  (see  his  equation  (96)),  and  yet  another  contribution 
due  to  ions  (his  equation  (97)).  Our  fluid  model  yields  the 
compressional  energy, 


2 

bw  - 


16n 


2 

25? 


Y+  V 


Zq  V^(0)  , 


(102) 


eff 


where  is  the  bounce-averaged  effective  collision  finequeacy. 
Note  that  when  =  0,  (102)  reduces  to  (97),  as  it  should.  (If 

»  y  •  Oien  the  factor  j/iy  +  can  be  approximated  by 
Y/Vjff).  Thus  b'^W g  (our  equation  (102))  includes  the 
compressional  effect  due  to  both  electrons  and  ions,  and  that 
there  is  no  separate  ion  contribution  as  Galeev's  equation  (97) 
suggests.  For  large  (or  equivalently,  large  stochastic 

diffusion),  the  stabilizing  effect  of  b^W^  can  be  strongly 
reduced.  Under  these  conditions,  it  is  possible,  in  principle,  to 
recover  an  electron  tearing  insubility,  but  there  is  no  ion  tearing. 
This  conclusion  supports  the  recent  results  of  Pellat  et  al.  [1991] 
who  question  the  very  existence  of  ion  tearing,  but  contradicts 
the  fmdings  of  Kuznetsova  and  Zelenyi  [1991]. 

In  view  of  the  controversy  in  analytical  theories,  much  can  be 
learned  from  particle  simulations.  Unfommately,  electromagnetic 
particle  simulations  of  collisionless  tearing  inevitably  involve 
making  compromising  choices  on  such  parameters  as  mjmi ,  the 
system  size  (which  determines  the  range  of  unsuble 
wavelengths),  and  the  spatial  grid  size.  We  have  cited  several 
such  simulations  earlier,  and  it  is  fair  to  say  that  in  all  of  them,  an 
insubility  with  the  theoretically  predicud  growth  rau  and 
characteristics  of  ion  tearing  has  been  very  difficult  to  futd. 
Since  we  believe  that  both  electron  and  ion  dynamics  (which  are 
tied  by  the  constraint  of  quasi-neutrality)  should  be  retained  in 
simulations  of  collisionless  tearing,  we  first  comment  on  reported 
results  from  two-species  simulations  that  include  a  5,  Geld. 
Swift  and  Allen  [1987,  P.10,013]  report  that  their  previous 
unpublished  work  showed  “no  evidence  of  the  development  of 
any  type  of  iiuubility.”  They  also  attribute  correctly  the 
observed  subility  to  ion  compressibility.  Zwingmann  el  al. 
[1990]  report  results  mostly  for  the  mass  ratio  mJmi  -  1  ,  with 
some  discussion  of  a  case  with  =  1/10  .  As  they  note,  the 

case  mJmi  =  1  caiutot  distinguish  between  electron  artd  ion 
tearing.  (If  an  ion  uaring  mode  exists,  its  growth  rau  should  be 
much  larger  than  the  electron  uaring  growth  rau  when  the  mass 
ratio  is  realistic.)  Their  resulu  show  significant  discrepancies 
with  theoretical  predictions  [Schindler,  1974].  In  particular,  the 
growth  rau  observed  in  the  simulation  is  up  to  an  order  of 
magnitude  less  than  predicud  by  theory.  We  attribuu  the 
growth  of  the  insubility  in  these  simulations  for  small  values  of 
5,  to  electron  uaring,  not  ion  Uaring.  This  hypothesis  can  be 
tested,  of  course,  by  a  study  which  compuus  the  growth  rau  as  a 
furution  of  mjin^  . 

Apart  from  two-species  simulations,  there  are  one-species 
simulations  of  the  ion  uaring  mode  in  which  the  electrons  tre 
involved  only  as  a  static  charge-neutralizing  background 
[Terasawa,  1981;  Hamilton  and  Eastwood,  1982;  Swift,  1983; 
Ambrosiano  and  Lee.  1983;  Pritchett  et  al.,  1991].  It  is  clear 
from  our  previous  discussion  that  these  simulations  cannot 
realistically  simulau  electron  tearing  modes.  Furthermore,  any 


inference  on  (he  viability  of  (he  ion  tearing  mode  from  these 
simulations  is  questionable  because  electron  dynamics  has  been 
eliminated  arbitrarily  for  reasons  of  computational  convenience. 

We  conclude  this  section  with  the  remark  that  unless  the  B, 
field  is  very  small,  field  lines  cannot  reconnect  to  form  islands  in 
the  linear  regime.  The  striking  contrast  between  configurations 
with  £,  =  0  which  tear  easily  to  form  magnetic  islands  and 
configurations  in  which  significant  values  of  S,  inhibit  tearing 
is  illustrated  well  by  Figures  6.2.4  and  6.2.9  in  Galeev’s  review 
paper.  In  Figure  6.2.4,  islands  develop  at  the  separatrix  where 
collisionless  reconnection  provides  accessibility  to  a  stale  of 
lower  energy.  In  Figure  6.1.9,  there  is  no  well-defined  seperatrix, 
and  the  system  sustains  global  compressional  oscillations. 

5.  Effects  of  akd  Spatially  Varying  fl, 

WBL  considered  the  effect  of  a  constant  By  field 
superimposed  on  the  two-dimensional  configuration  of  section  4. 
Their  analysis  of  the  electron  tearing  mode  dealt  with  the  inner 
region  dynamics  using  kinetic  theory,  but  global  aspects  of  the 
dynamics  such  as  the  bounce  motion  was  neglected.  The  aim  of 
the  present  effort  is  to  explore  the  consequences  of  these  global 
effects  in  the  context  of  the  improved  asymptotic  equilibria 
developed  in  section  2. 

At  first,  we  consider  equilibria  in  which  By  is  constant,  and 
the  spatial  dependencies  of  B^  and  B^  are  described  in  region 
1  (near-Earth)  by  (18)  and  (19).  in  region  2  (middle)  by  (20)  and 
(21),  and  in  region  3  (distant-tail)  by  (22)  and  (23).  respectively. 
As  noted  in  section  2.  the  spatial  structure  of  B^  in  the  near- 
Earth  and  middle  regions  are  similar,  except  that  the  magnetic 
field  is  weaker  and  (he  current  sheet  is  wider  in  the  middle 
region.  It  is  shown  in  Appendix  A  that  the  bounce  period  in 
region  2  is  larger  by  an  order  of  magnitude  than  in  region  1. 
In  region  3.  since  fi,  is  vanishingly  small,  Xj,  is  extremely 
large. 

Certam  conditions  must  be  fulfilled  for  the  electron  tearing 
instability  to  occur.  First,  there  must  be  magnetic  free  energy 
available  to  drive  the  insubility;  i.e.,  we  must  have  <  0  , 
where  is  given  by  (89).  This  means  that  the  stability 

parameter  A,'  must  be  positive.  In  standard  analyses  of 
collisionless  subility  of  the  tail  (see,  for  instance,  GaUev  [1984, 
and  references  therein]).  A,'  is  replaced  by  Ag' .  Note  (hat  this 
overestimates  the  range  of  unstable  wavelengths  because  A2' 
(equation  (78))  is  negative. 

Second,  the  subilking  compressional  energy  5^^  due  to  the 
bounce  triotion  of  electrons  should  not  exceed  the  desubilizing 
term  .  We  show,  a  posteriori,  that  in  region  1,  yr^  «  I  , 
but  in  regions  2  and  3,  we  have  2  1  and  yxj,  »  1  , 
respectively.  The  compressional  stabilization  is  thus  significant 
in  region  1  but  not  so  in  regions  2  and  3. 

In  regions  2  and  3,  the  effect  of  the  electron  bounce  can  be 
neglected.  We  are  then  back  in  the  framework  of  WBL  who 
obtained  the  dispersion  equation  for  electron  tearing  modes 
neglecting  electron  bounce.  From  a  kinetic  analysis,  carried  out 
in  Appendix  B.  we  obtain  the  complex  frequency  to  =  tOj^  iy , 
where 
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in  region  2.  Here  o).  =  ifc  ,  jq  is  <  constant  which  can  be 
shown  to  have  a  numerical  value  of  approximately  10  under 
typical  conditions,  L,  =  and  A/=  ieJikAg 

X^)  «  a/  .  For  comparison,  we  recall  that  if  =  0 ,  the  real  and 
imaginary  parts  of  the  complex  growth  rate  are,  respectively 
[Drake  and  Lee,  1977],  given  by 
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In  (103)  and  (104),  all  equilibrium  quantities  are  evaluated  at  r  = 
0,  In  particular,  since  v(0)  /  £(0)  =  BylB^iz  =  0).  we  can  write 
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where  a  n  By/Bg  .  If  a  -  £  ,  a  necessary  condition  for 
insubility  is 

a;  =  Ag'  +  Aj'  >  0  .  (108) 

which  gives 

itX<l/V2.  (109) 

Equation  (109)  implies  that  wavelengths  larger  than  2V2KX  C= 
9  Rg  for  X=  1/lf)  may  be  unstable.  In  order  for  the  instability 
to  grow,  however,  it  must  also  satisfy  the  condition  yr^  ^  1 .  A 
viable  class  of  instabilities  is  obtained  for  2t  3X  ;  these  do 
obey  the  condition  >  1  for  sg  =  10 ,  =  2X/s  . 

If  By  «  B, ,  i.e.,  £  »  a  ,  the  condition  A/  >  0  for  instability 
reduces  to  XX  <  a/£  .  As  shown  by  WBL,  this  condition 
predicts  unstable  wavelengths  which  are  much  too  large  to 
account  for  reconnection  events  in  the  near-Earth  and  middle 
regions. 

We  note  that  (104)  and  (106)  have  been  obtained  from  a 
kinetic  analysis,  and  except  for  a  factor  of  ,  can  also  be 
obtained  from  the  fluid  model  The  fluid  analog  of  (106)  is  (74), 
derived  in  section  3.  The  fluid  analog  of  (104)  has  been  derived 
in  Appendix  B,  and  is  given  by 

Jh'.  S.  ’ 

Y=  A,.  (110) 

2*0  2., 

As  before,  the  results  from  the  fluid  and  kinetic  calculations  differ 
by  a  multiplicative  factor  of 

In  regions  where  yr^  ^  1  ,  the  desubilizing  effect  of  By  may 
be  understood  as  follows.  If  By  >  0 ,  Caleev  [1984]  poinu  out 
that  the  energy  spent  for  plasma  compression  is  the  work  done 


by  the  perturbed  plasma  current  generated  by  the  perturbed 
pressure  gradient. 
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Now,  if  1,  then  the  bounce  average  of  (112)  should  be 

taken,  and  the  second  term  on  the  right-hand  side  of  (112) 
averages  to  zero.  On  the  other  hand,  if  yx^  21 .  or  »  1, 
then  there  is  no  bounce  average  to  be  taken,  and  the  second  term 
tends  to  reduce  the  first  term  because  T\jl  J -  Bj IBt . 

There  is  some  evidence  in  the  numerical  simulations  of  Swift 
and  Allen  [1987]  that  the  presence  of  By  enhances  the  tearing 
activity  near  z  s  0  compared  with  the  case  By-0  .  (See  their 
section  4J.)  Clearly,  there  is  a  need  for  two-species  simulations 
including  By  using  either  asymptotic  equilibria  of  the  kind 
developed  in  this  paper,  or  numerical  solutions  of  the  equilibrium 
Grad-Shafranov  equation  [cf.  Voigt  and  Hilmer,  1987], 

Finally,  we  comment  on  equilibria  with  spatially  varying  By, 
discussed  in  section  2.2.  There  we  show  that,  for  a  class  of 
pressure  profiles,  the  spatial  variation  of  By  can  cause  the 
formadon  of  an  X  point  on  the  r  =  0  line.  We  show,  furthermore, 
that  this  can  occur  at  near-Earth  distances  for  average  values  of 
By  and  B^  characterisdc  of  the  plasma  sheet  The  configuration 
thus  formed  is  likely  to  be  highly  magnetically  stressed.  Under 
such  condidons,  rapid  reconnecdon  may  occur  at  the  separatrix 
in  both  the  linear  and  nonlinear  regimes.  We  conjecture  that  the 
collisionless  reconnecdon  rate  in  this  geometry  is  likely  to  be 
much  larger  than  the  rates  derived  in  this  paper.  Such  a  geometry 
calls  for  a  separate  treatment  and  the  exploration  of  that 
possibility  is  left  to  future  work. 


6.  Conclusions 


This  paper  makes  two  main  contribudons  to  the  problem  of 
collisionless  tearing  modes  in  the  Earth’s  magnetotail.  The  first 
involves  the  development  of  asymptodc  magnetotail  equilibria 
including  all  three  components  of  the  magnetic  Geld,  with  realism 
in  the  modeling  of  the  normal  component  of  the  magnetic  Geld,  G, 
(x,z) .  The  second  involves  the  development  of  a  fluid  model  that 
is  physically  transparent  and  accurate  in  reproducing  the 
parametric  dependencies  of  the  growth  rates  of  collisionless 
tearing  modes  calculated  from  kinetic  theory. 

One  of  the  significant  conclusions  of  this  paper  is  that  the  ion 
tearing  mode,  which  has  been  the  subject  of  considerable  research 
and  controversy  over  the  last  two  decades,  does  not  occur.  This 
is  true  for  both  two-  and  three-component  models  of  the 
magnetotail.  We  are  not  the  first  to  suggest  this,  because 
Lembege  and  Pellat  [1982]  and  Pellat  et  al.  [1991]  have  preceded 
us,  albeit  in  the  context  of  the  simple  two-component  equilibrium 
(1).  We  find  that  for  the  two  and  three-component  equilibria 
given  in  this  paper  that  if  there  is  a  collisionless  tearing 
instability  in  ihe  magneutail.  it  is  the  electron  tearing  mode. 


There  are  cenain  condiiions  that  must  be  fulfilled  for  the 
electron  tearing  mode  to  be  seen.  A  significant  value  of  A,  in 
the  two-component  magnetotail,  represented  by  (1),  tends  to 
suppress  the  instability.  The  reason  for  this  strong  subilizaiion 
can  be  understood  in  dynamical  terms.  In  our  view,  the  dynamics 
are  a  symptom  of  a  deeper  cause  which  has  to  do  with  geometry. 
The  main  cause  of  the  subilizadon  of  the  tearing  mode  is  that 
destroys  the  separatrix  at  r  =  0  .  By  contrast,  if  we  set  =  0 
but  include  a  Bj  field,  the  separatrix  at  z  «  0  is  undisturbed, 
and  election  tearing  modes  can  easily  occur. 

The  main  difficulty  posed  by  the  magnetotail  is  that  all  three 
components  of  the  magnetic  field  can  be  signiCcaiu.  Then,  the 
circumstances  that  favor  electron  tearing  are  those  that  minimize 
global  features  of  the  dynamics  such  as  electron  bounce,  and  keep 
the  electron  confined  near  z  s  0  .  It  is  intuidvety  clear  that  the 
By  field  tends  to  confine  electrons  near  z  s  0  .  and  hence  helps 
the  electron  tearing  instability  grow.  The  asymptotic  equilibria 
presented  in  this  paper  have  regions  where  the  stabilizing  effect 
of  electron  bounce  can  be  neglected  and  where  the  presence  of  By 
-  B,  can  cause  the  excitation  of  electron  tearing.  A  tearing 
iiuubility  in  which  field  lines  actually  undergo  genoine 
topological  change  does  not  occur  unless  B«  is  very  small. 
Unless  topological  change  occurs,  the  instability  is  likely  to 
saturate  nonlinearly  at  a  relatively  low  amplitude.  We  do  not 
believe  that  such  a  weak  instability  can  account  for  the  dramatic 
signatures  associated  with  current  disrupdon  and  diversion 
during  substoims. 

The  insubility  is  more  interesting  when  B,  is  zero.  We  have 
demonstrated  that  if  we  include  By  and  allow  it  to  vary  spatially 
in  a  three-dimensional  magnetotail  equilibrium,  then  B,can 
vanish  at  near-Eanh  distances.  The  linear  as  well  as  the  nonlinear 
growth  of  electron  tearing  modes  in  such  a  configuradon  is  likely 
to  lead  to  interesting  results  and  will  be  invesdgated  in  the  near 
future. 

An  important  challenge  for  a  theory  of  substoims  is  that  it 
should  account  not  only  for  the  violent  acdvity  that  is  associated 
with  substorms,  but  also  idendfy  conditions  under  which  the 
magnetotail  is  stable.  A  universal  instability  that  occurs  always 
and  spontaneously  is  likely  not  to  be  a  correct  explanation 
because  that  would  suggest  the  magnetotail  is  always  unstable, 
which  is  not  observed  to  be  the  case.  In  this  work,  we  have 
identified  conditions  under  which  electron  tearing  modes  may  be 
unstable  and  delineated  regimes  when  they  are  noL  It  is  our  hope 
that  this  paper,  as  well  as  its  forerunner  (WBL),  will  stimulate  a 
reexaminadon  of  old  as  well  as  new  data  in  substoims  with  a 
renewed  emphasis  on  the  By  field.  Observadonalists,  many  of 
whom  we  have  cited  here,  have  been  aware  over  the  last  13  years 
of  the  ubiquitous  presence  of  the  By  field,  varying  spatially  as 
well  as  in  time  before  and  during  different  phases  of  a  subsiorm. 
What  is  required  is  a  more  systemadc  study  coirelating  By  and 
B,  with  the  occurrence  of  substorm  onset.  More  two-species 
electromagnetic  particle  simulations,  including  all  three 
components  of  the  magnedc  Geld,  are  also  required,  both  to  check 
analytic  theory  and  to  model  realistically  global  features  of 
magnetotail  equilibria. 


APPENDIX  A:  Bounce  Motion  Of  Electrons 


3 


An  electron  gyrating  along  a  field  line  in  the  magnetouil  may. 
under  certain  conditions,  bounce  between  two  turning  points. 
Since  (x,  z)  »  B„(x,  -  z) ,  the  z  coordinates  of  the  turning 
points  may  be  written  as  z  =  ±  z, ,  where  the  constant  z,  is 
determined  by  the  parallel  energy  of  the  electron. 
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where  p  is  the  magnetic  moment  The  bounce  period  is  defined  as 


where  /  is  the  coordinate  along  a  field  line.  Using  the  relation 
dl  I B  =  dUB I  ,  we  have 


I, 

J  w  if,  (1-1^, 


We  recall  that  5,=  tanh  z  and  fl,  =  e(l  -  z  tanh  z) .  For  I  z I  £ 
1  ,  we  use  the  approximations  fl,  =  z  and  fl,  r  £(1  -  z*)  . 
Since  ^  1)  =  ^  •  consider  only  the  domain  0  ^ 

I  z  I  <  1 .  Equation  (A3)  may  be  then  approximated  as 

f'  r  2  2  ,  2.2, 1/2 
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Defining  u2  =  I  -  z/z, .  we  get 


If  we  take  B^=t  everywhere,  as  in  Lembege  and  Pellat  [1982], 
we  get 

^  {  2  16  2 
^4  =  —  J(l-«)d«  =  -S— z^  (A6) 
6  0  3c  v.^ 

For  z,  -  1  ,  equation  ( A6)  predicts  a  bounce  period  of  =  22  s 
(for  c  -  0.1  and  v„  -  2/s) .  We  show  below  that  this  is  much 
smaller  than  the  bounce  period  in  our  asymptotic  equilibria. 

The  integral  (AS)  can  be  evaluated  exactly  to  give 
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Clearly,  as  I  z,l  t ;  this  mearu  that  some  electrons  are 

lost  and  some  have  very  large  bounce  periods.  An  average 
bounce  period  for  the  confined  electrons  can  be  obuined  by 
averaging  over  the  distribution  function  of  electrons.  For  fixed 
z  we  have 


i 


where  is  the  maximum  parallel  velocity  above  which  the 
electron  is  lost.  Note  that  (A8)  underestimates  the  bounce  period 
because  it  should  be  normalized  by  the  fraction  of  confined 
elections  which  is  smaller  than  1  .  Since  all  bouncing  electrons 
pass  through  z  s  0,  we  take  the  distribution  at  z  =  0  to  be 
Maxwellian.  We  estimate  ■  v,,  =  l)/e  -  1]*^  . 

Asymptotic  evaluation  of  the  integral  (A8)  gives  <  ‘r^>  > 
3««/ev^.  For  e  r  0.1,  v„  =  2/z  ,  we  get  <  t*  >  •■  10*j  in 
region  1.  In  region  2,  since  £«  is  space-dependent,  e  is  replaced 
by  e'  ■  e~  e/4  ,  and  the  length  scale  is  amplified  by  (see 
section  (2.1)).  Then  <%>  ia  amplified  by  a  numerical  factor  of 
approximately  which  gives  <Ts>  2  lO^z .  Fory^  KT^x  , 
we  thus  have  y<x^>  «  1  in  region  1,  but  y<x^>  i.  1  in 
region  2. 


AppendocB: 

Dispersion  Relation  For  Electron  Tearing  Modem 

In  order  to  keep  this  paper  self-contained,  we  review  here  the 
derivation  of  the  dispersion  relation  for  electron  tearing  modes 
using  kinetic  theory.  The  main  effort  lies  in  calculating  the 
perturbed  current  from  the  perturbed  distribution  function 
/la  by  means  of  the  relation 


It  is  convenient  to  use  (32)  to  write 


fu 


We  recall  that  [Laval  a  ai,  196d;  Galeev,  1984] 


r’  f  ®/a  c  1 


2. 


Hence,  from  Ampere’s  law, 

J,  =  V^A,  =  So  I 

1  4jt  1  n 


we  get 


L  iL^ 

dx^*dz^ 


X,  cosh  (z/X.) 


^l'*'  ,2 


\  cosh  (zJX) 
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dy  v  f 
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From  the  linearized  Vlasov  equation  it  follows  that 


-7 

dl  ■'la 


E.-v/ 


(Bl) 


(B2) 


(B3) 


(B4) 


Oa' 


(B5) 


(B6) 


where  dtdt  m  d/dt  -t-  v  •  V  and  /og,  is  the  Maxwellian 
distribution.  Integrating  (B6)  along  characteristics  and 
substituting  the  result  in  (B5X  we  get 

fil  ilh  2  A 

xWrz/X) 


-7  Z  I?-  jdy  y(t)j  y(tyE^(nfojn  dt' 


Alt 

m  - —  o  •  E,  ,  (B7) 

c  1 

where  ^  is  the  collisionless  conductivity  tensor  which  can  also 
be  written  as  (see,  for  instance,  Horton  and  Tajma  [1990]), 


dt  y(0  y(t  - 1) 


^  «*P  J  • 


-ik  j  \ 


Near  r  =  0  .  in  the  inner  region  l9/3zl  »  1 3/3x1  -  lA .  Then 
a  can  be  diagonalized  in  the  form 


o  0  0 


a  =  0  a  0 


0  0  o,. 


where  x  =  x.  t=a,y  -a^r,  bn  B/S  =  a^y  +  a,  t  with 
and  a,=  B,lB.  We  v^te  £ix  ^  ot,£i^ -a^i,. 
Ell  =  ayEiy  + a,Ei,,  since  a  +a,  =  1. 

Defining  V]  m  Ai^ ,  Xi  ■  ^ii  •  wd  writing,  as  in  WBL,  any 
perturbed  quantity  such  as  V  in  the  form  \|/(z)  exp(ikx  -  iat)  , 
we  get 

c  A  “ 

—  =  y  •  o  •  E, 

4ji  ♦-* 


.  (ajoj  +  a* op  £  +  otyOt.fo^  -  op  £,  ,  (BIO) 
where  Ox  ■  Oj »  Og .  Similarly, 

C  A  /  2  2  \Z 


dz  /  2  2  \r 


(BID 


The  parallel  conductivicy  0]  can  be  calculated  by  using  the  drift- 
lunetic  approximation  for  elecffons.  It  can  be  shown  that  [Drake 
and  Ue.  1977;  WBL  1990] 

2 

J,  =  i  —H  (0)  -  (0*)  /  r(s)  .  (B 12) 

“  4Jt  <B  ' 

where 

2 

(0  2 

o.  =  i — Z(s)  .  (B13) 

“  4*0) 

where  tn*  =  JkV^, ,  j  =  (o/(i:||  y^) ,  «nd  Z(r)  is  the  plasma 
dispersion  function. 

The  perpendicular  conductivity  is  mostly  due  to  the  ion 

polarization  drift,  since  the  £  x  B  drift  carries  no  current.  The 
full  ion  response  [Cowley  et  al.  1986;  WBL]  essentially  reduces 
to 

-  nmr^  - 

J,  =— i-(-io)E,)  =  -i(0 - tE,.  (B14) 

■^5  47iV: 


where  =  Bti^nn  mi)^^  is  the  Alfv6n  velocity.  Note  that 
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for  magnetotail  plasmas.  We  then  recover  the  result  derived  by 
WBU  i.e.. 
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(B16) 


We  now  introduce  two  characteristic  frequencies.  One  is  (a«a 
.  The  other  is  a>e  a  where  is  the  time  it  takes  for 
an  electron  to  travel  a  characteristic  distance  Ir^  along  x 
[WBL].  From  the  field-line  equation 


dz 

B 

f 


(B17) 


it  follows  that  if  an  electron  travels  along  x,  it  must  travel 
along  z ,  where  z^  is  given  by  the  relation 

z5  =  X^l-exp(2£/a)].  (B18) 


We  then  obtain  [WBL] 


34 


T 


=  txp  (eJkX).  (B19) 

2c  v„ 

We  incroduce  the  strained  inner  variable 

4  =  ^  a  ^  .  (B20) 

tf 

where  5«  >  ((H/k  y,,)  and  f o  =  1  ^  to^  >  (Oe  (as  in 

region  3),  but  sq  =  if  o)^  <  a>c  (as  in  region  2).  The  shear 
length  L,  is  defined  by  L,m\  B(z  »  0)  Bq  .  With  these 
definidoRS,  (B16)  can  be  rewritten  as 
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where  =  Cos5,(0)5(j.  and 

(1  -  CO*/©)  5]  kl  , 

H  -  Z  (§'). 

which  is  the  same  as  (77)  of  WBL, 

Equation  (B21)  can  be  diagonalized  as 


where 


=  “oV  V’ 

Oj  -  -e^v  +  Ogi. 


(B21) 
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(B23) 

(B24) 

(B25) 


As  shown  by  WBL,  02  (X  ®  0)  “  0  •  ’•'Wch  yields 

¥(0)  =  (aycg)  x(0).  (B26) 

The  soludon  for  a\  is  different  in  different  regioiu. 

In  region  1,  where  ‘yt^«  1 .  the  integradon  over  z  involves  a 
bounce  average  which  extends  over  the  outer  region.  We  can 
then  see  the  subilizing  effect  of  this  bounce  by  noting  that  it 
tends  to  cancel  out  the  perturbed  current  in  the  inner  region  and 
gives  ¥  =  i  *  0 .  In  regions  2  and  3,  however,  we  have  yt*  >  1 
and  yTj  »  1  ,  respecdvely.  Integrating  (B23)  over  the  inner 


region,  and  using  the  "constant  v"  approximation,  we  get  the 
jump  condition 

^  =  (00^  +  4)  0,(0)  (B27) 

““  ia.«r 

In  order  to  match  the  inner  region  to  the  outer  region,  we  use  the 
relations 

'  d  ~  1 0*  d  -  I+- 

A  «  j:  log  v  L.  =  77  log  V  . .  (B28) 

1  dZ  outer  '  ^  dt  outer  ' 

*  A  mm 

A  *s  T"  log  X  Ll  =  '7~  log  Z  (B29) 

2  dz  *^oo«r'0-  dz  **am«tr'-- 
The  left-hand  side  of  (B27)  can  then  be  written  as 

da.  .  da.  da.  , 

'm«r  *  d^.  •  d^  '  — 

aocr  ^  outer 

=  5^  [og  VfO)  A,  +  eg  x(0)  Aj  ]  .  (B30) 

The  right-hand  side  of  (B27)  gives 
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Using  (B30)  and  (B31).  we  obtain  the  dispersion  equation 

£>(t.(0)  +  fl-— lz(0)  =  0.  (B32) 

where 


A' 


2  ' 
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2  2 


(B33) 


Equation  (B32)  is  slightly  different  from  the  analogous  equation 
derived  by  WBL.  The  difference  can  be  attrributed  primarily  to 
a  slightly  different  treatment  of  the  asymptotic  matching 
condition:  here  we  have  extended  the  outer  solution  all  the  way  to 
z  B  0  instead  of  taking  it  to  the  outer  limia  of  the  inner  region. 
In  the  rotation  of  WBL.  A'-Ar  +  iA/;  however,  if  we  recall  that 
A/  -  e  «  A^ .  it  is  clear  that  the  results  derived  by  WBL  and  here 
do  not  differ  significantly.  (In  order  to  avoid  possible 
confusion,  we  point  out  that  <0*  has  been  neglected  in  certain 
results  derived  in  WBL;  this  affects  the  real  part  of  to  .  but  not 
the  growth  rate.)  Solving  the  dispersion  equation  (B32).  we 
obtain  (O^  and  y  ,  given,  respectively,  by  (103)  and  (104)  in 
region  2  and  (lOS)  and  (106)  in  region  3. 

The  growth  rate  derived  from  kinetic  theory  can  also  be 
obtained  from  the  fluid  model.  The  fluid  analog  of  (106)  has 


already  been  derived  in  section  3.  Here  we  derive  easily  the 
anolog  of  (104).  Recalling  the  heuristic  discussion  given  by 
WBL  (see  section  3.  WBL),  have  y-  (Oga  = 

(n^kdgKXoi,).  Here  d,  =  AV2kQ  .  which  is  obtained  from  (68). 
with  given  by  (88).  If  we  now  use  (B25).  we  get  (110) 
which  is  larger  than  the  growth  rate  derived  from  kinetic  theory 
by  the  factor  . 
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